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Abstract 

We investigate the scenario where the dark matter only interacts with the charged leptons in the 
standard model via a neutral vector mediator Z'. Such a scenario with a 430 GeV dark matter can fit 
the recent positron fluxes observed by the AMS-02 Collaborations, with the reasonable boost factors. 
We study the possibility of searching such leptophilic Z' via its lepton final states and invisible decay 
modes at the future electron-positron colliders, such as the International Linear Collider (ILC) and 
the Compact Linear Collider (CLIC). We find that for the benchmark models with Z' mass from 
1.0 TeV to 1.5 TeV, the searches for the invisible decays of Z' —> \X is easily achieved at the CLIC 
1.5 TeV runs via the mono-photon process. However, lighter Z' with mass from 0.5 TeV to 0.8 TeV 
are challenging to see. The di-lepton plus single photon channel can reveal the Z' mass at the ILC 
and CLIC with moderate luminosities. 
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I. INTRODUCTION 


The existence of dark matter (DM) has been confirmed by many cosmological and as- 
trophysical observations, while little is known about the nature of the DM candidates and 
their interactions with the Standard Model (SM) particles in the visible sector. Experimen¬ 
tal searches to the DM candidates and their interactions with the SM particles include the 
underground direct searches, indirect searches and the collider searches. Among them, the 
most intriguing results are the recent observations of the excessive high-energy positrons that 
cannot be explained by the known astrophysical sources. One possible explanation for these 
excessive high-energy positrons is due to the DM contributions through annihilation or decay 
processes. The excess of the high-energy positron fluxes was first observed by PAMELA m- 
[3]. The most recent measurement of the positron energy spectrum extends the positron 
kinematic energy up to 300 GeV [3j. Later, the observation of the excessive high-energy 
positrons was reported by the Alpha Magnetic Spectrometer (AMS) collaboration [1], More 
recently, the AMS-02 extended the measurements of the positron fraction up to ~ 500 GeV 
with improved precision M- It is also shown that the positron fraction no longer exhibits 
an increase with energy above 200 GeV [7J- 

There have been a lot of efforts to address the anomalous positron fluxes observed by 
the PAMELA and AMS-02 results, including the Sommerfeld enhancement effects iKsi, 
the lightest neutralino as DM candidate in the supersymmetric models mm, and other 
possibilities [2TL433j . Here we consider a leptophilic DM setup with a Z' vector boson being 
the mediator between the visible sector and the dark sector. In this setup, the DM candidate 
annihilates into three generations of charged leptons equally through the Z' interactions. 
We do not consider the DM annihilation channels to quarks or gluons since there are no 
excessive anti-proton observed in the PAMELA and AMS-02 yet [SSI 135| . The current direct 
search experiments and hadron collider searches also provide constraints to the DM effective 
couplings with quarks and gluons [36HS0] , With a boost factor of BF ~ 0(100), the DM 
candidate with mass of 430 GeV is likely to fit the AMS-02 observation of the excessive 
positron fluxes at the high energy region. This leptophilic Z' together with its coupling 
to the DM candidate can be tested at the future high-energy e + e _ colliders, such as the 
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international linear collider (ILC) [3T| and the compact linear collider (CLIC) [32]. The 
direct discovery of such leptophilic Z' can be most easily achieved through the opposite-sign- 
same-flavor (OSSF) di-lepton final state associated with a photon. We analyze the potential 
to discover the interaction between Z' and DM in this setup at the future CLIC runs at 
1.5 TeV. Our results show that for Mz< (1-0,1.5) TeV, the mono-photon channel appears 
quite promising to detect the interaction after imposing a set of appropriate kinematic cuts. 

This paper is organized as follows. In section [D} we describe the setup of the simplified 
leptophilic DM model, where a vector boson Z' mediates the DM interactions with the charged 
leptons via the s—channel. A set of benchmark models are listed where the Z' vector boson 
masses are chosen to be Mz' G (0.5 TeV, 1.5 TeV). We also discuss the relic densities for 
the leptophilic DM. In section HI we discuss the indirect detections of DM candidates via 
the positron fluxes in the leptophilic Z' model, and compare with the latest AMS-02 results. 


Section IV is devoted to analyze the discovery potential of the leptophilic Z' at the future 
e + e~ colliders such as ILC and CLIC. The discovery of the leptophilic Z' is easily achieved 
via the OSSF di-lepton plus single photon final state, with the integrated luminosities of 
f Cdt ~ 0(1) pb -1 — 0(10) fb _1 for the benchmark models. Furthermore, we analyze the 
searches for the DM at the CLIC 1.5 TeV runs via the mono-photon channel. Finally, we 
conclude in section EO 


II. THE MODEL SETUP 

A. The leptophilic Z' with the DM couplings 

For models with extended gauge symmetries such as U(l) and/or SU(2), a neutral vector 
boson Z l is quite usual as a new physical state. Assuming a DM particle in the spectrum na¬ 
ns], it is natural to consider the Z' as the s —channel mediator between the visible and the 
dark sector. Such a Z l is likely to couple with different DM candidates with various spin 
natures. The details of the Z l couplings with different DM candidates and their annihilation 
cross sections can be found in Ref. |46j. For our illustration in this work, we assume a Dirac 
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fermion x as the DM candidate with the general couplings with the Z' as follows, 


£dm = (9x + 9*75 )x z 'p ■ 


( 1 ) 


Several previous studies have carried analysis of the direct, indirect, and the collider searches 
for the Z l interactions with the DM candidates dMZ]. In Refs. pHSaET], the mediator Z l 
is assumed to couple with both SM quarks and leptons. Consequently, the constraints from 
the direct searches for DM should apply. There also exist strong constraints from the LHC 
searches for the mediator via the dijet process [58] and the DM candidates via the mono-jet 
process m respectively. To avoid such constraints from mono-jet process and to address 
the possible positron excesses from the latest AMS-02 results, we assume a leptophilic Z' 
coupling of the following form 


u = Y , 1 + 9 / 75 )^;, ^ = 1 , 2 , 3 , ( 2 ) 

i 

with li representing three generations of charged leptons. In the sequential SM case for the 
Z' vector boson, the vectorial and axial couplings in Eq. ([ 2 ]) are given by 


9i = 


-(1 -4s 


4c w v ~ ~ wJ ’ ~4c w ’ 

with the short notation of the Weinberg angle (sw , cw) = (sin Q\y , cos 0\v)- In our discussions 
below, we simplify the Z l couplings with the DM and the charged leptons as follows 


9e = 


9 


(3) 


A x = 9x = 9x i x e = 9 X = 9 / , (4) 

in Eqs. 0 and 0. The A x and \i characterize the overall Z' coupling strengths with the 
DM candidate % and the charged leptons respectively. 


B. The experimental constraints 


There are several experimental constraints to consider for the leptophilic Z' couplings given 
in Eq. ([ 2 ]). The leptophilic Z' couplings will contribute to the lepton anomalous magnetic 
dipole EH] as 


A(9 - 2)e ~ 


6 vr 2 M|, ' 


(5) 
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Specifically, the measurement of the muon anomalous magnetic moment from the E821 ex¬ 
periment at BNL m shows a discrepancy j62] with the SM predictions of 


Aa^ = a® xp - a® M = 288(63)(49) x 10~ n , (6) 

where = ( g M — 2)/2. Accordingly, the upper bound of the coupling A^ from the g — 2 
constraint is 

^<5.6x10-3(1^). (7) 

This constraint is moderate for the very heavy Z' with mass of Alz' ~ 0(1) TeV. 

The coupling of \i is mostly constrained by the LEP-II searches for the Z'. For a vectorial 
coupling to electrons, this bound reads j63j 

A e < 0.25 x (^) , (8a) 

\AMe < 0.14 x (4^) for l + e. (8b) 

The Z' will be considered in the mass range of Aiz' G (0.5 TeV, 1.5 TeV) below. According 
to the LEP-II constraint listed in Eqs. ([8]), we fix the Z' coupling with charged leptons as: 

f 0.05, Aiz' < 1.0 TeV, 

Xe={ (9) 

0.1, Alz' > 1.0 TeV, 

and vary the Z' coupling with the DM A x only. 


C. The DM relic density 


The relic densities of the DM candidates are obtained by solving the Boltzmann equation 
for the evolution of the DM number density n 

— = -3 Hn--(av)(n 2 -n 2 eq ), (10) 

where n eq denotes the DM equilibrium number density, H is the Hubble parameter, v is the 
“relative velocity”, and a is the usual spin-averaged cross section 1 . It is well-known that the 


1 The factor of \ in Eq. (101 is due to the Dirac fermion as DM. 
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thermally averaged annihilation cross section (av) plays the key role in determining the DM 


relic densities. The solution of Eq. (10) yields the DM relic density as 

2 2.14 x 10 9 GeV’ 1 

VLh z ~ — 


Jgl ,2 M pl 


( 11 ) 


where h is today’s Hubble parameter, M p \ = 1.22 x 10 19 GeV, and g * is the number of 
effective degrees of freedom at the freeze-out temperature. The quantity J is given by 

° (<™> 


J = 


-dx. 


( 12 ) 


where x = m/T and Xf is the value given at the freeze-out temperature. In the non-relativistic 
limit, the thermally averaged annihilation cross section is given by 

2x 3 / 2 


(av) ~ — j=r- I (av)e 1 ^ 2 exp(— xe)de , 

/?r Jo 


(13) 


with e = (s — 4m 2 )/(4m 2 ) and 


1 


av = 


647r 2 (s — 2m^.) 


1 - 


4m 


/ 


dQ \M \ 2 , 


(14) 


where \M\ 2 is the spin-averaged amplitude square. Far from the resonance and particle 
production threshold, one can often expand the annihilation cross sections in powers of e, 


av « a + be + 0(e 2 ) 


Making use of Eq. (13), one obtains the thermally averaged annihilation cross section 


(av) ~ a+-bx 1 + 0(x ). 


(15) 


(16) 


Thus, when the annihilation occurs not close to a resonance, the relic density can be approx¬ 
imated by 


2 2.14 x 10 9 GeV -1 x f 

Qh ~ -— . 

gl /2 M p i(a + 3b/4xf) 


(17) 


However, when the annihilation is near the resonance, one should determine J(xf) by per¬ 
forming the numerical integral of 

r°° v 2 (av\ 

J(xf) = / dv —-— / dxx^ 1 / 2 exp(— xv 2 /4). (18) 

Jo v 47r J Xf 
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The details of the near resonance DM annihilations are given in Refs. [641 [65 ]. 

In our setup, there are two primary DM annihilation channels depending on the mass 
relations between Mz> and m x : 

1. For m x < Mz'-, the DM annihilate into the SM lepton pairs through the s—channel Z' 
exchanges. 

2. For m x > Mz>, the DM either annihilate into the SM lepton pairs through the 
s— channel Z' exchanges; or annihilate into the Z' pairs through the t— channels. 


M z , [ TeV] 

A* 

Dcdm^ 2 

BF 

r vis ( GeV) 

r inv ( GeV) 

0.5 

1.1 

0.12 

154 

0.2 

0 

0.6 

0.8 

0.12 

175 

0.2 

0 

0.7 

0.6 

0.11 

134 

0.3 

0 

0.8 

0.24 

0.12 

134 

0.3 

0 

1.0 

0.20 

0.13 

330 

1.6 

0.8 

1.1 

0.44 

0.12 

223 

1.8 

6.0 

1.2 

0.68 

0.12 

207 

1.9 

17.7 

1.3 

0.90 

0.13 

217 

2.1 

37.3 

1.4 

1.23 

0.12 

191 

2.2 

79.7 

1.5 

1.52 

0.12 

192 

2.4 

138.4 


TABLE I: The benchmark models with Dirac DM mass of m x = 430 GeV and different M z > inputs. 
The couplings of A x are chosen so that the corresponding relic densities are Dc DM /i 2 ~ 0.12. The BF 
denotes the boost factor needed to fit the positron fluxes of the AMS-02 experiments. The last two 


columns are the decay widths of Z', calculated by Eqs. (28) and (30). 


In our analysis, we only consider the cases where the DM annihilations occur not close 
to the resonance. Hence, the approximate results of Eq. is sufficient for evaluating the 
DM relic densities. To get the benchmark models, one should consider the DM relic densities 
which are related to the DM annihilation cross sections. Making use of the approximate 
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(19) 


expression of the relic density in Eq. (|17[), we have 


Qh 2 ~ 


1.07 x 10 9 GeV -1 x f tt(M|, - 4 m 2 x ) 2 
y/giM p i l + 5/xf 3A|A |m| 


where a phase space factor S reads 

1 + 20 ml/Ml, , 

S ~ 4(1 -4 /M|,) ' (20) 

The current measurements of the cold DM relic density from both WMAP and Planck 
satellite m give Dcdm^ 2 — 0.12, which will be considered as the major constraint for us 
to obtain the benchmark models. In Table. |T| we list the leptophilic DM benchmark models 
for the further studies. Here, we choose the DM mass as m x = 430 GeV in order to fit the 
current AMS-02 positron flux spectrum, to be discussed in the next section. The masses 
of Z' vector boson are taken in the range of Myj e (0.5 TeV,1.5 TeV), where we avoid 
the near-resonance choice of My' = 0.9 TeV. The coupling strength A x between Z' and 
DM for each case is taken as in Table. [I] so that the corresponding DM relic density fits the 
observed value of Qh 2 ~ 0.12. The evaluation of the DM relic densities here are made by 
micrOMEGAs m- The details of other elements in the table will be given in the following 
context of the indirect detections and the collider searches. 


III. THE INDIRECT DETECTIONS OF DM 

In this section, we discuss the indirect searches for the DM candidates following the model 
setup above. In our case, the DM candidate annihilates into three generations of charged 
lepton pairs equally. Besides of the e + e~ final states, the (/r + /x“ , t + t ~) final states also 
contribute to the positron energy spectrum due to their decays. Some details can be found in 
Refs. [7l,[72]. The more recent work on fitting the AMS-02 results into various leptonic final 
states include Refs. [731476] . where the DM annihilations into different leptonic final states 
were considered separately. The positron propagation in the Milky Way is described by the 
following diffusion equation m- 

V • [K(E , r)Vf+\ + [b(E , r)/+] + Q(E , f) = 0 , 


( 21 ) 




where /+ is the number density of positron per unit energy, K(E , r) is the diffusion coefficient, 
b(E , r) is the rate of the energy loss, and Q(E , r) describes the source of positron. The 
primary positron sources are due to the DM annihilation 


Q(E,r) 


®’5>> 


dN + 

~dE 


( 22 ) 


where the summation goes through all possible channels that generate positrons in the final 
states. The primary positron flux from the DM annihilations is given by 

>*T(E) = ^f+(E,r e ), (23) 

where r® ~ 8.5 kpc is the distance between the Sun and the galactic center. In practice, we 
obtain the CalcHEP [68] model files for the benchmark model by using the LanHEP package [69] . 
These model hies will be passed to micrOMEGAs [70] in order to evaluate the primary positron 
fluxes. The DM density p(r) is chosen to be the local density at the Sun p® « 0.3 GeV/cm 3 
and the halo prohle function F(r ): 


P(r) = Pe F ( r ) 


(24) 


where the default Zhao prohle m in the mi cr OMEGAs is taken. In addition, there are sec¬ 
ondary positrons due to interactions between cosmic rays and interstellar medium, whose 
hux can be well approximated as EZmiHD]: 

4.5 E 0 - 7 


&%°(E) = - 


+ 650 E 2 - 3 + 1500 E 4 - 2 


GeV • crri • sec • sr 


-i 


(25) 


With the primary and secondary positron huxes given in Eqs. (23) and (25) respectively, 
we combine them into the total positron hux as 

&f( E ) = $+ ime {E) x BF + $^ C (F), (26) 

where a boost factor (BF) is included. This BF parametrizes the effects of local DM clumps 


distributions in the area where DM annihilation takes place. As shown in Eq. (22), enhancing 
the local density by a factor of 10 would yield an enhancement of the positron hux by a factor 
of 100. Afterwards, we perform the following y 2 ht 


X 2 = 


E 


( p” odel - /+ AMS ) 2 


(T? 


(27) 
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FIG. 1: The predictions of positron fluxes from the \X —> Z' —> £ + £ processes for different DM 
masses and the M z > = 1 TeV case together with the AMS-02 measurements of the positron fluxes [5]. 


to the AMS-02 data of positron fluxes in order to determine the corresponding BF for each 
benchmark model in Table. [T| The summation in Eq. (27) goes through each energy bin for 
the positron spectrum from 10 GeV up to 430 GeV. Below ~ 10 GeV, the counting rate of 
the positrons becomes unstable, which is affected by the solar modulation [6]. Here /xf odel 


are the products of E 3 • with the total positron fluxes given in Eq. (26). /x^ MS and 

cji are the AMS-02 measurement of the positron fluxes E 3 • 4* and errors in the energy range 
of E £ (10,430) GeV. As an example, for the m x = 430 GeV and Mz> = 1 TeV, we get the 
minimum of y 2 /d.o.f ~ 0.87 when BF = 330. We demonstrate the positron flux predictions 
from the benchmark model with m x = 430 GeV and Mz> = 1 TeV in Fig. [lj The models 
with different values of Mz> predict the same positron flux once some appropriate BF’s are 
taken into account. To fit the AMS-02 predictions, we expect a large boost factor (BF ~ 100) 
for large mass of DM. The fits to other benchmark models and the corresponding BF’s are 
listed in Table. [T] previously. 

Now, we give some comments on our fit to the AMS-02 positron flux. First, our fit is 
based on the assumption that the Z l couples with the three generations of leptons equally. In 
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general, the positron spectrums generated from the e, //, r final states have different shapes. 
If the fraction of each lepton in the final state is not fixed, but determined after fitting to 
the data, then there would be a much better agreement between the experimental data and 
the theoretical prediction from DM annihilation [MS]. Second, we notice that the present 
experimental result only cover the energy spectrum of the positron flux up to about 500 GeV. 
As a result, the DM of a mass of about 430 GeV would be able to fit this spectrum. At the 
same time, because of the mass limit, the positron spectrum induced by the DM is ended at 
about 430 GeV. In the future, it is likely that the AMS-02 data show a continuing existence 
of the positron flux over a larger energy range. In order to account for this possibility, we 
also show the fit with DM masses of 600 and 800 GeV in Fig. [I] It is seen that they start to 
deviate from the present data, which means that the DM of a mass heavier than 600 GeV 
with equal couplings to leptons is less likely to explain the current AMS-02 data in this 
leptophilic Z' model. 

IV. THE COLLIDER SEARCHES FOR Z' AND ITS INVISIBLE DECAYS 

The leptophilic DM with Z l interaction is just one of the many models trying to explain 
the positron fluxes observed. This explanation can be tested by other experiments, such as 
the collider searches. Given that the vector boson Z' is not too heavy, it is likely to produce 
Z l directly at the LHC and the next generation e + e _ colliders. The search for leptophilic Z' 
at the hadron colliders is challenging due to the large QCD background, where the jets can 
be misidentified as leptons. In this work, we focus on the search at the next generation e + e~ 
colliders. We expect that not only a heavy Z' is discovered, but also its coupling with the 
DM can be confirmed. 
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A. The decays of Z' and its searches at the e + e colliders 


The vector boson Z' couples to both SM leptons and Dirac DM The partial decay 
width to the SM leptons is 


T[Z'^£ + £-] = 


Mz' 

~Y2tt 


' 1 - 


4 rrij 

Mf, 




m 2 , - 1 


'Ml 


M 2 Z , 


(28) 


Under the massless limit of me = 0 and the simplified couplings of A t = gj = gf, one has 
this visible partial decay widths become 


T [Z’^tT] = ^\j. 

07r 


(29) 


The invisible decay width to the Dirac DM x reads: 





(^) 2 (l + 2^f) + (^) 2 (l- 4mX 


Ml 'J 


2 , 
) 


M'b 


(30) 


with the coupling simplifications of A x = = 9\ i n the seconc l One. In Fig. [ 2 J we show 

the decay branching fractions BR[Z 7 ] by combining the visible and invisible decay widths in 
Eqs. (28) and (30). The evaluations are made for the benchmark models given in Table. [Tj 
For the mass range of Mz> E (0.5 TeV, 1.0 TeV), the decay mode of Z' —> £ + £~ dominates. 
When the invisible decay channel of Z' —> XX is kinematically allowed, i.e., Mz> > 2m x , the 


invisible decay mode becomes more dominant. 

Below, we will study the leptophilic Z' searches at the future high-energy e + e - colliders 
based on the benchmark models that we listed in the previous Table. [T| Considering the 
mass ranges and interaction properties of Z 7 , we are interested in the e + e _ colliders with 
center-of-mass energy yfs > 1 TeV. Therefore, we explore the experimental searches at 
the ILC 1 TeV run and the CLIC 1.5 TeV run. In practice, we implement the benchmark 
models into UFO model hies by using FeynRules [81] for the event generation. Events are 
generated for both signals and the corresponding SM background processes by Madgraph [52]. 
Afterwards, the PYTHIA [83] and PGS [84] will be used for the parton shower, hadronization, 
and fast detector simulations. For the detector simulations, we implement the design of 
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FIG. 2: The BR \Z' —► t + t ] (red curve) and BR \Z' —> \x] (blue curve) for Mz> G 

(0.5 TeV, 1.5 TeV). 

the ILD [85] in the PGS. For both ILC and CLIC runs, we consider the polarization scheme 
of P(e + ,e~) = (—30%,+80%) for the incoming electron and positron beams, which will 
suppress the SM background contributions efficiently [ 86 ], [87] 

B. The direct searches for the Z' through the di-lepton plus photon channel 

We first discuss the searches for Z' through the OSSF di-lepton associated with a single 
photon channel: e + e _ —> Z\-^r £ + £~)^, where t = (e,/i). In particular, we require that 
the Z' is produced on-shell. For the benchmark models we listed in Table. [TJ we will study 
the signals for the M 7 j G (0.5 TeV, 0.8 TeV) samples at the ILC 1 TeV run, and the 
Mz' G (1.0 TeV, 1.5 TeV) samples at the CLIC 1.5 TeV runs. The relevant SM background 
processes contributing to the OSSF di-lepton plus photon final states include: e + e~ —> £ + £~ 7 , 
e + e _ —> IT + LF _ 7 , and e + e _ — > r + T~ r j. It turns out only the first process of e + e _ —> £ + £~7 
contribute dominantly, while the cross sections for all other processes are sub-dominant. In 
addition, events containing jets and large Ifir from the e + e _ — > W + W~ 7 and e + e~ — > t + t~ 7 
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processes can be vetoed, hence the contributions from these two processes can be safely 
neglected. 

For illustration, we list a set of kinematic cuts to be imposed at the ILC 1 TeV runs for 
the benchmark models with Mz' € (0.5 TeV, 0.8 TeV). 

• Cut-1: all events containing jets with px > 10 GeV and \r]\ < 2.5 are vetoed. Besides, 
we also veto events with fix > 50 GeV. 

• Cut-2: the events containing OSSF di-lepton plus a photon are selected. The OSSF 
di-leptons should satisfy the requirements: |r/(t' ± )| < 2.5, px{?o) > 20 GeV, px{h) > 
10 GeV, and 2.0 < A R(£ + £~) < 5.0. Here and l\ represent the leading and sub¬ 
leading charged leptons ordered by their transverse momenta. 


• Cut-3: The photon in the selected events satisfy the requirements: |??( 7 )| < 2.5 and 
pr(n) > 20 GeV. Since the photon is produced via the e + e - —> yZ' process, the 
photon energy is determined by (s — M|,)/2- v /s. Considering the detector smearing 
effects, we impose the following cuts to the single photon energy E^: 

S — My, S — My, 

z - 25 GeV < E 1 < z + 25 GeV . (31) 

2ys 2ys 

The corresponding E-y distributions for the signal and the SM background processes 
are shown in the left-panel of Fig. [3} 


• Cut-4: finally we reconstruct the invariant mass of the OSSF di-leptons around the 
Z' mass window: \Mpp — Mz>\ < 50 GeV. The distributions of the invariant mass of 
OSSF di-leptons Mpp for the signal and the SM background processes are shown in the 
right-panel of Fig. [3j 


To search for the benchmark models with Mz> E (1.0 TeV, 1.5 TeV) at the CLIC 1.5 TeV 
runs, we shall follow the similar cut flows listed above for the ILC 1 TeV runs. Meanwhile, 


we modify the single photon energy cut window from 25 GeV in Eq. (31) to 50 GeV. Fur¬ 
thermore, the cut to the invariant mass window of OSSF di-leptons Mpp will be modified to 
|Mpp — Mz’\ < 100 GeV for the CLIC 1.5 TeV runs. 
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ILC at ts=1.0 TeV 


ILC at ts=1.0 TeV 




FIG. 3: The distributions of the single photon energy E 7 (left panel) and the OSSF di-lepton invariant 
masses (right panel) for the signal process of e + e“ —► Z'{—> l + £~)7 and the SM background process 
of e + e“ —> £ + £~"/ at the ILC 1 TeV runs. 



cut -0 

cut -1 

cut -2 

cut-3 

cut-4 

5 G 

M Z ' = 0.5 TeV 

138 

124 

118 

114 

111 

145 pb " 1 

e+e-j 

1,125 

1,050 

619 

39 

31 

- 

M z , = 0.6 TeV 

215 

191 

187 

180 

175 

65 pb -1 

+ 

1 

1,125 

1,050 

619 

40 

35 

- 

M Z ' = 0.7 TeV 

238 

210 

205 

198 

192 

61 pb -1 

+ 

1 

1,125 

1,050 

619 

46 

41 

- 

M z , = 0.8 TeV 

447 

395 

384 

373 

360 

27 pb -1 

£+f"7 

1,125 

1,050 

619 

59 

54 

- 


TABLE II: The cross sections (unit: fb) after imposing kinematic cuts for the benchmark models 
from the OSSF £ + t~ plus single photon channel and the relevant SM background processes at the 
ILC 1.0 TeV runs. In the last column for different benchmark model signals, we list the integrated 
luminosity needed for the 5 <r excess via the OSSF di-lepton channel. 


After imposing the series of kinematic cuts, we list the cut efficiencies for the benchmark 
models at the ILC 1 TeV runs and at the CLIC 1.5 TeV runs in Table. |TT] and Table. IIIII 
respectively. The integrated luminosities required to observe 5er excess from the OSSF di¬ 
lepton plus a single photon channel for each benchmark model are also listed. For benchmark 
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cut-0 

cut-1 

cut-2 

cut-3 

cut-4 

5 a 

M Z ' = 1.0 TeV 

265 

228 

222 

216 

211 

43 pb" 1 

£+l~ 7 

594 

553 

339 

29 

27 

- 

M Z ' = 1.1 TeV 

119 

101 

99 

93 

91 

159 pb -1 

0V”7 

594 

553 

313 

30 

28 

- 

M z , = 1.2 TeV 

78 

65 

63 

57 

55 

409 pb -1 

£+£~ 7 

594 

553 

313 

35 

33 

- 

M z > = 1.3 TeV 

69 

57 

55 

46 

44 

762 pb -1 

+ 

1 

594 

553 

313 

48 

45 

- 

M Z ' = 1.4 TeV 

66 

54 

51 

37 

35 

1.55 fb” 1 

£ + £~j 

594 

553 

313 

70 

64 

- 

M z , = 1.5 TeV 

25 

20 

19 

4 

4 

22.5 fb” 1 

£ + £~ 7 

594 

553 

313 

15 

13 

- 


TABLE III: The cross sections (unit: fb) after imposing kinematic cuts for the benchmark models 
from the OSSF £ + £~ plus single photon channel and the relevant SM background processes at the 
CLIC 1.5 TeV runs. In the last column for different benchmark model signals, we list the integrated 
luminosity needed for the 5 a excess via the OSSF di-lepton channel. 

models of Mz> £ (0.5 TeV,0.8 TeV), the Z l can be searched via the e + e” —> Z\-^r l + £~)7 
channel when the 1 TeV ILC runs with the integrated luminosities reaching 0(10) pb” 1 . For 
benchmark models of Mz> £ (1.0 TeV , 1.5 TeV), the signal channel of e + e” —> Z\-^r £ + £~) r ) 
can reach 5(7 C.L. at the 1.5 TeV CLIC runs with the integrated luminosities reaching 
0(10) pb” 1 -0(10) fb" 1 . 

C. The direct searches for the invisible modes of Z mono-photon channel 

Besides the direct discovery of the leptophilic Z' at the future e + e” colliders, here we 
discuss the search for the invisible decay modes of Z' via the mono-photon process of 
e + e” — > Z'{—} xx)l- The previous studies on the DM searches at the e + e” colliders and the 
implications to the DM direct and indirect searches were carried out in Refs. |86l 1881(91] . For 
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CLIC atVs=1.5TeV 


CLIC atVs=1.5 TeV 




FIG. 4: The distributions of the fp (left-panel) and the single photon energy E 1 (right-panel) in the 
mono-photon signal process of e + e _ —► Z'(—> xx)l and the corresponding SM background processes. 
The benchmark model with Mz> = 500 GeV is taken for the CLIC 1.5 TeV runs. 


the mono-photon channel, we should look for the 7 + pr signal. The corresponding SM back¬ 
ground should be: e + e _ — > 7 + 7 , which are due to the e + e~ — > Z(—t vv)x and the t— channel 
IT 1 * 1 —exchanging processes. Other processes like e + e _ — y £ + £~ r y and e + e _ — y t + t~ 7 are also 
likely to fake the 7 + tfsp signal, when neither of the charged leptons or jets are tagged in the 
detector. The cross sections for the relevant SM background processes at the CLIC 1.5 TeV 
runs with the beam polarizations of P(e + , e~) = (—30% , +80%) read 


-7 ui"y] = 408 fb , 

(32a) 

-» £ + £~ 7 ] = 594 fb , 

(32b) 

—> t + t~ 7 ] = 13 fb . 

(32c) 


Again, the cross sections are evaluated by considering the preliminary cuts imposed at the 
parton level by Madgraph. 

Now we discuss imposing the appropriate kinematic cuts in order to select signal events 
from the SM background processes. For the benchmark models with 500 GeV < Mz> < 
800 GeV, the mono-photon process of e + e~ —> \X1 is produced through the off-shell ex¬ 
changes of the Z’ . Accordingly, we illustrate the set of kinematic cuts to be imposed: 

• Cut-1: the events containing jets and leptons with px > 10 GeV and \r]\ < 2.5 are 
vetoed. The events containing 7 + lf}x signal will be selected, where the single photon 
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cut-0 

cut-1 

cut-2 

cut-3 

5(7 

M z > =0.5 TeV 

9.7 

9.4 

9.3 

9.3 

113 fb” 1 

M z , = 0.6 TeV 

6.1 

5.9 

5.8 

5.8 

290 fb” 1 

M z > =0.7 TeV 

4.3 

4.2 

4.1 

4.1 

580 fb” 1 

M z > = 0.8 TeV 

0.99 

0.98 

0.94 

0.94 

11 ab” 1 

z/^7 

408 

404 

401 

390 

- 


TABLE IV: The cross sections (unit: fb) after each kinematic cut to be imposed for the benchmark 
models in the mono-photon process and the relevant SM background processes at the CLIC 1.5 TeV 


runs. 


should further satisfy the conditions: pt(,i) > 10 GeV and ( 77 ( 7 ) | < 2.5. 


• Cut-2: the events with Ifix < 400 GeV are selected, and the corresponding distributions 
of the Ifrr for both signals and SM background processes are shown in the left panel of 
Fig-0 after imposing the Cut-1. 


• Cut-3: the single photon energies in this set of benchmark models should be bounded 
from above: < (s — 4m^,)/(2- v /s) ~ 503 GeV for the CLIC 1.5 TeV runs. The 

corresponding distributions of the single photon energy E * y for both signals and SM 
background processes are displayed in the right panel of Fig.0after imposing the Cut-1 
and Cut-2. We also note an equivalent kinematic variable to consider is the recoil mass 
of: 


TTlrec 



pj) 2 - 


(33) 


For the set of benchmark models with M z > E (0.5 TeV, 0.8 TeV), one expects that 
m rec > 2 m x = 860 GeV. The cut of the recoil mass turns out to be equivalent to the 
Ery cut here. 


Afterwards, we list the kinematic cut efficiencies in Table. [TV| after imposing Cut-1 to Cut-3 
for the benchmark models with Mz' E (0.5 TeV , 0.8 TeV). The searches of the mono-photon 
channel for this set of benchmark models require the integrated luminosities of ~ 0 ( 100 ) fb ” 1 
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or even larger, which are generally challenging. This is due to the smallness of the mono- 
photon cross section with an off-shell invisible decay mode Z' —>• XX produced. 



cut-0 

cut-1 

cut-2 

cut-3 

cut-4 

5(7 

M z > = 1.0 TeV 

240.6 

240.0 

239.9 

233.5 

231.6 

21 pb^ 1 

Z/Z/7 

408 

404 

402 

9.3 

6.1 

- 

M Z ' = 1.1 TeV 

672 

665 

664 

619 

608 

7 pb _1 

PZ/7 

408 

404 

401 

14.1 

10.3 

- 

M z > = 1.2 TeV 

1,173 

1,153 

1,148 

1,022 

1,001 

4 pb” 1 

PZ/7 

408 

404 

397 

23.8 

19.0 

- 

M z > = 1.3 TeV 

2,010 

1,963 

1,926 

1,643 

1,608 

3 pb^ 1 

Z/Z/7 

408 

404 

387 

48.8 

42.2 

- 

M Z ' = 1.4 TeV 

3,819 

3,703 

3,681 

3,202 

3,147 

2 pb” 1 

Z/Z/7 

408 

404 

375 

161 

150 

- 

M z , = 1.5 TeV 

2,264 

2,108 

2,084 

1,759 

1,747 

4 pb -1 

Pi/7 

408 

404 

348 

161 

158 

- 


TABLE V: The cross sections (unit: fb) after each kinematic cut to be imposed for the benchmark 
models of M Z ' £ (1.0 TeV,1.5 TeV) in the mono-photon process and the relevant SM background 
processes at the CLIC 1.5 TeV runs. 


For the benchmark models with Mz> £ (1.0 TeV, 1.5 TeV), the invisible decays of Z' —> 
XX occurs on-shell. We impose the following set of cuts to distinguish the signals from the 
SM background: 

• Cut-1: the events with 7 + pr are selected. The selected photon should satisfy the 
conditions: 177(7) | < 2.5, pr{l) > 10 GeV. The events with leptons or jets of > 

10 GeV and \p(£,j)\ < 2.5 are vetoed. 
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Cut-2: we set an upper bound to the pr of the selected events: 


M Z ' - 

= 1.0 

TeV : 

: 0 

GeV 

< 

Pt 

< 

450 

GeV, 

(34a) 

M Z ' - 

= 1.1 

TeV : 

: 0 

GeV 

< 

Pr 

< 

400 

GeV, 

(34b) 

M Z ' - 

= 1.2 

TeV : 

: 0 

GeV 

< 

Pr 

< 

300 

GeV, 

(34c) 

M Z i = 

= 1.3 

TeV : 

: 0 

GeV 

< 

Pr 

< 

200 

GeV, 

(34d) 

M Z t - 

= 1.4 

TeV : 

: 0 

GeV 

< 

Pr 

< 

150 

GeV, 

(34e) 

M/i - 

= 1.5 

TeV : 

: 0 

GeV 

< 

Pt 

< 

100 

GeV. 

(34f) 


• Cut-3: the single photon energy in the signal events is still determined by E 7 = 
(s — M|,)/2i/s. As for the CLIC 1.5 TeV case, we determine the single photon energy 
as: E™ e (E 7 - 25 GeV , E 7 + 25 GeV). 


Cut-4: the recoil mass m rec defined in (33) can be used to reconstruct the Z' mass. 
Thus, we require the cut of |m rec — M z < | < 50 GeV for the M Z i E (1 TeV, 1.5 TeV) 
samples. 


In Table. [VJ we list the cross sections for the 7 + pr signals and the dominant SM background 
of vv'y after each kinematic cut imposed. For the CLIC running at \Js = 1.5 TeV, our analysis 
shows that the mono-photon signal searches for the benchmark models can reach 5 a excess 
with the integrated luminosities of 0(1) pb” 1 . 


V. CONCLUSIONS 

The current indirect DM search results are intriguing for probing the DM natures and 
their possible connection to the visible sector described by the SM. The search for the cos¬ 
mic positrons from both PAMELA and the recent AMS-02 collaborations keep showing the 
disagreement between the measurements and the known astrophysical predictions. The DM 
candidates annihilating into various SM final states were considered as possible source for 
such excessive high-energy positron fluxes. In this work, we hypothesize the couplings be¬ 
tween Dirac DM particle x and the charged leptons in the SM mediated by a leptophilic 
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Z l 2 3 4 5 . After considering the experimental constraints from muon anomalous magnetic moment, 
and the LEP-II searches, we show that this setup is capable of fitting the current AMS-02 
results of the positron fluxes. 

Such a setup with a leptophilic Z' can be clearly tested in the future high-energy e + e~ 
collider, with the center-of-mass energy running at 0(1) TeV. The searches for the Z' is 
considered for the OSSF di-lepton plus a single photon final state. We analyze the signals 
for the benchmark models and the corresponding SM background contributions, and find 
the discovery of the Z' requires moderate luminosities at ILC/CLIC runs. Furthermore, the 
searches for the DM productions are essential for determining the DM couplings to Z l . For the 
benchmark models with Z' mass of Mz r G (1.0 TeV, 1.5 TeV), the searches for the invisible 
decays of Z' -A XX is easily achieved at the CLIC 1.5 TeV runs via the mono-photon process. 
Meanwhile, for the light mediator cases of Mz> G (0.5 TeV , 0.8 TeV), the mono-photon cross 
sections are typically small, making the discovery of the invisible mode challenging. 
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